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cavities, while the current rigorous quantum treatment and numerical assessment sheds 
light on the feasibility of experiments that aim at optically manipulating beams of 
massive molecules with multimode cavities. 
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1. Non-resonant, cavity-assisted laser cooling 

Laser cooling of atoms has been an important driving force in quantum optics 
throughout the last three decades [H [21 [3]. It has led to great advances in atomic 
spectroscopy and atomic clocks [4] , atom interferometry [5] , and quantum sensors P, [7] , 
and it even paved the way towards the new physics of quantum degenerate gases [HI [9] . 
The optical manipulation of atoms has become a routine tool in many experiments 
world-wide. A limitation of laser cooling arises, however, when we need to address 
systems with a more complex internal structure than atoms. Large molecules, atomic 
clusters or solids in general do not exhibit suitable dipole-allowed transitions that could 
be used for running a closed cycle of absorption and spontaneous reemission to the 
initial state, as required for the controlled dissipation of kinetic energy into a photonic 
emission. Laser cooling of particle motion based on spontaneous light forces is therefore 
mostly precluded. 

It has, however, been realized that even then optical manipulation and cooling 
can be achieved [10] when it is done on particles in a high-finesse optical cavity that 
preferentially redirects any incident scattered radiation because of the enhanced mode 
density in a given direction. Blue-shifted Rayleigh scattering of light at a particle that is 
moving along the cavity axis may thus open a channel for damping the particle's kinetic 
energy. The cooling rate then depends primarily on the real part of the particle's off- 
resonant polarizability x well as on the light intensity and its frequency detuning 
with respect to the cavity resonance. In a standard setting the incident laser light is 
red-detuned {ojl < ujc) with respect to a cavity resonance such that the Doppler shift 
of a particle moving towards one of the cavity mirrors will support the preferential 
emission into the cavity mode. The combination of a blue-shifted photon recoil with 
a velocity-dependent enhancement of the scattering rate leads to a friction force that 
effectively decelerates the particle. 

Resonator-assisted optical cooling has been theoretically proposed for the strong 
coupling regime [10] and extended towards the collective cooling of a cloud of 
particles [HI [121 [131 [El [IS] o^' the cooling of complex molecules, too [TH [17] . 
The schemes were extended from the one-dimensional case with different pumping 
configurations [HI [191 [20] to the three-dimensional case, still using a single cavity [2T] . 
Even a lowering of the internal temperature of the particles was studied in a cavity- 
mediated environment |22j. 

As of today, the idea of non-resonant cavity-assisted optical cooling has been 
experimentally demonstrated for individual atoms [23l [24] as well as for cold atomic 
ensembles. Collective effects in the interaction between an atomic cloud and the cavity 
field led to the discovery of collective recoil lasing [251 126] . 

A related field has been opened recently in a complementary setting: The laser- 
assisted cooling of micromechanical oscillators [271 [2H] shares much of the physics and 
the formalism with cavity- assisted cooling of atoms [29j. The manipulation of cantilevers 
may thus be understood as a specific form of a more general scheme in which one may 
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either study the free motion of individual atoms or the bending of a micromechanical 
cantilever that is restricted to a harmonic oscillation close to the quantum mechanical 
ground state. As tailor-made micromechanical oscillators couple to the light very 
efficiently and as they withstand high laser intensities, the cooling powers in such 
experiments were shown to be substantial (see [30], and references therein). 

A number of experiments were already successful at both ends of the mass scale, 
either with atoms or using small solids. It is, however, still an open challenge to 
apply cavity cooling schemes to freely moving mesoscopic objects. Such nanoparticles 
could for instance be the Ceo-fullerenes [31], heavier biomolecules [32] or large 
silica nanospheres [33l |3l] which have recently been used or proposed for quantum 
superposition experiments. The implementation of efficient cooling and slowing methods 
will be a key issue in all those experiments. 

On the other hand, most of these examples will operate with objects that have 
broad optical lines and low oscillator strengths, without being large enough to form a 
mirror for visible or even infrared light. All these cases have therefore to be regarded in 
the weak-coupling limit [35], as outlined below. 

In this article we present a full quantum description of the motion of such a 
massive polarizable particle under the influence of a strong laser field and a possibly 
large manifold of high-finesse cavity modes. This reproduces known results from former 
semiclassical derivations but also exceeds those works, as we provide experimental case 
studies for the motional damping of various nanocluster particles and a proper treatment 
of the strong pump enhancement of the coupling to the cavity modes. We furthermore 
show how the presence of the multimode cavity leads to a velocity damping effect 
acting on the particle, which increases with the number of accessible high-finesse modes. 
An explicit calculation of this effect is carried out on the basis of a realistic confocal 
cavity setup with an accessible manifold of ~ 10^ degenerate modes. We predict an 
enhancement of the optical friction force by more than two orders of magnitude when 
compared to a single-mode cavity. The experimental parameters are assessed using 
clusters with masses up to 10^ amu, where we also include the absorption and the 
free-space scattering of photons into our model. 

In the first section we start by deriving a master equation for the motion of a 
polarizable particle in a pumped cavity mode. The explicit derivation is carried out 
with a single particle in a single cavity mode. We then generalize the model to a one- 
dimensional treatment of particles in M modes. In Section |3] the master equation 
is translated into phase space in order to identify the effective friction force and the 
diffusion terms that act on the particle to lowest order in the photon recoil. An explicit 
formula for the velocity dependence of the force beyond linear friction is presented 
there. We incorporate absorption and Rayleigh scattering into our model in Section H] 
and discuss its validity with respect to experimental realizations with large particles, 
such as atomic clusters. The results are assessed in a parameter regime which is realistic 
for those mesoscopic particles. 
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2. Master equation for a particle in a pumped cavity 

We start by considering a polarizable particle placed in a single-mode high-finesse cavity 
which is characterized by its resonance frequency Uc and the photon loss rate 2k. A 
laser at frequency ujp pumps the resonator continuously. 

The coupling between the particle and the cavity field is mediated through the 
ofi^-resonant scalar electric polarizability x ^ which we assume to be real, positive and 
constant within the frequency range between Up and Uc - and the quadrature components 
of the quantized electric field operator 



We denote operators by non-serif letters throughout this article. Here, is the mode 
volume and f{r) the mode function fulfilling the Helmholtz equation Af{r) = —k'^f{r) 
with the wave number k = 2n/X. For a cavity with widely extended plane-parallel 
mirrors the mode function is simply f{x) = cos /ex. For more practical finite cavities 
and curved mirrors the transverse beam profile is closer to Gaussian and the mode 
function is three-dimensional, including both the mode profile and phase shifts related 
to the wavefront curvature [36] . 

The interaction between the particle and the cavity is described by the field energy 
of the induced dipole in rot at ing-wave- approximation. 



where we have introduced the coupling constant Uq = —UcXf^^oVc- It has the dimension 
of a frequency and describes both the optical potential created by a single photon and 
the shift of the cavity resonance due to the presence of the particle, i.e. the change of 
index of refraction. Note that if the particle motion is one- dimensional and confined 
around a position where the mode function can be linearized, |/(a;)p — )■ x, one recovers 
the well known optomechanical interaction Hamiltonian [271 [28]. 

We describe the interaction in a rotating frame with respect to the pump frequency 
Up and in a displaced frame with regard to the stationary pump amplitude a such that 
the driving Hamiltonian for the empty cavity vanishes [271 EH]- The mode operator 
transforms according to a — > a + a. As opposed to previous treatments based on weak 
field assumptions [35], we assume the presence of a strong pump, i.e. a significant 
population of the pumping mode |a| ^ 1 and introduce a first weak coupling condition: 

The coupling to the field should be small compared to the cavity damping rate, 
|f/o| ^ so that the presence of the particle introduces only small fiuctuations of the 
cavity field around the value a. 

We may therefore neglect the term in the displaced coupling Hamiltonian that is 
quadratic in the displaced mode operators, and we approximate the total cavity-particle 
Hamiltonian with 




(1) 




(2) 




p 



(3) 
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where A = — cUp is the frequency detuning between the pump and the cavity. 
The Hamiltonian is the sum of the free cavity energy He = /iAa^a, the energy 
of the particle that is moving in the optical potential generated by the stationary 
pump Hp = p^/2mp + /i?7o|a/(r)p, and a linear coupling to the field fluctuations 
H/ = /if/o|/(r)pQ;*a + h.c. Together with the cavity damping this yields the master 
equation for the density operator p describing the motion of the polarizable particle and 
the cavity field 

dtp = ~ [H, p] + K (2apat - atap - pa^a) = (£p + + Cc) p. (4) 

It is written in terms of the Liouville superoperators Cpjp = —i[Hpj,p]/h, and 

Ccp= — zA[a''"a, p] + /t(2apa''" — a"''ap — pa''"a). (5) 

Our goal is now to derive an effective master equation of the particle motion by 
eliminating the field degrees of freedom. In order to be allowed to adiabatically eliminate 
the cavity we require the second, stronger weak coupling condition: 

The effective linear coupling to the field fluctuations is assumed to be small 
compared to the cavity damping rate, \Uoa\ <^ k. This means that any excitation 
of the cavity on top of the pump amplitude a that is added by a scattering process at 
the particle will rapidly decay, and (a'^^a) -C 1 in the displaced frame. As the effective 
coupling strength is now enhanced by the pump field this represents a much stricter 
requirement than the first weak coupling condition, and it can be violated in practice 
by sufficiently strong pump fields, as we will see later. 

The second weak coupling condition allows us to use the projection formalism 
[37t [38[ l39] to eliminate the cavity state and derive an effective master equation for the 
particle. Following this procedure, we introduce the orthogonal superprojector V onto 
the cavity vacuum state (relative to the stationary pump), Vp = trc(p) ® |0)(0|. It 
fulfills = V, and its complement is denoted by Q = id — P. The goal is to find 
a closed evolution equation for Vp, which yields the master equation for the reduced 
particle state pp = trc(p) after tracing out the cavity. Both projectors V, Q only 
act on the cavity subspace, and therefore commute with the particle's Liouvillian Cp. 
Furthermore, one may easily check that VCc = '^c'P = VCjV = 0. This allows us to 
split the master equation (jlj) into two coupled equations for Vp and Qp by applying the 
superprojectors from the left, 

dtVp = VCpVp + VCjQp, (6) 

dtQp= Q{Cp + Cc)Qp+Q'CiVp+QCiQp. (7) 

Here, the initial condition is an empty cavity at t = 0, Vp{0) = pp(0) ^ |0)(0| 
and Qp(0) = 0. The weak coupling condition implies that the contribution of field 
fluctuations to the state Qp and to the time evolution Cj are suppressed by the 
perturbation parameter \UQa\/K, ^ 1 on the fast decay time scale We solve the 

coupled equations up to first order in this parameter to obtain 

dtVp{t)=VCpVp{t)+ [ dTVCiQe''^^QCi{T)Vp{t), (8) 

Jo 
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with Cj{t) = e'^^'^/^/e"^^^ the r-delayed interaction LiouviUian. We approximate the 
upper bound of the integral by t ~ oo for t ^ 1/k. This is justified because of the 
fast field damping in the integrand related to the cavity evolution term e^*^"^. Equation 
([8]) has to be traced over the cavity to obtain an effective master equation for the 
reduced particle state pp{t). After insertion of all the Liouvillians, and several steps of 
straightforward calculation, the effective master equation reads 

dtPP = ~ [Hp,pp] - f/o>r 0/('-)l',gPP - PPgT (9) 



with the operator 



g 



In the following it will be convenient to refer to this object as the memory operator. It 
describes the delayed cavity reaction towards the past trajectory of the particle, 

|/(r,p;r)|^ = e-^H--/V(•■)|V^--/^ (11) 

within the memory time constant We will discuss below that the delay 

transformation can be approximated by the shearing |/(r, p;r)p ^ |/(r — pr/mp)|^, 
if the particle moves sufficiently fast and the influence of the optical potential of the 
pumped mode is negligible within I/k. Although we started our derivation from a 
fully Markovian treatment of particle and cavity field, the memory operator renders 
the effective master equation for the particle non-Lindbladian within the limits of the 
underlying weak coupling assumptions. This will later result in small negativities of 
the momentum diffusion coefficient for moving particles, which can be safely neglected 
in practice. We will furthermore see that it is the delayed response of the field on the 
particle motion that is responsible for the optical friction force - a phenomenon that has 
been well studied and experimentally observed [27l [28] with optomechanical systems. 
There the one-dimensional motion of a mechanical oscillator is confined at a position 
where the mode function, that mediates the field coupling, can be linearized. In the 
more general case of freely moving particles considered here the analoguous friction 
effect will depend on the particle trajectories through the cavity but can nevertheless 
be observed. Compare this e.g. to polarization gradient cooling [IQ], where it is the 
delayed response of the internal particle dynamics on the local field polarization which 
is responsible for the friction force. 

A generalization of master equation ([9]) to particles in M modes is given in 



Appendix A[ In the simultaneous presence of A^ particles in the same pumped mode. 



our weak coupling criterion reads, however, A^|f/oa| ^ k, and it may be violated by 
a beam or cloud of weakly coupling highly polarizable particles of sufficient density. 
Our criterion implies, in particular, that any cross-talk between the particles via field 
fluctuations in the cavity is small. It therefore also excludes collective correlation effects 
and self-organization flEl [26l [T3] . At the same time, we do not impose any restriction 
onto the number M of cavity modes. Using a cavity with many nearly-degenerate 
modes, such as in a confocal setup [HI HI], one may significantly enhance the control of 
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Figure 1. Scheme of a two-dimensional mode setup with one j/-mode strongly pumped 
to the stationary amplitude a, and standing- wave cavity modes orthogonally oriented 
to the pump along x. The latter are damped by the rate k and detuned from the 
pump. A particle with polarizability x experiences dissipative dynamics by scattering 
photons from the pump to the damped and detuned cavity. 



the particle motion. Furthermore, the generahzed result in Appendix A only requires 
one of the possibly many modes with similar damping and coupling strength to be 
strongly pumped by |a| ^ 1. Our general master equation therefore also includes two- 
dimensional setups, where the pump light illuminates the particles perpendicularly to 
the cavity axis [21]. The master equation of a single particle interacting with M empty 
modes, and one mode /o pumped at amplitude a, is given by 

M 

dtPP = — [Hp,pp] -J2U>\' ([/o(0/n('-),gnPp] + h.C.) . (12) 

n=0 

We here assume the polarizability x of the particle to be the same for all modes, i.e. 
Un = — X V'i^oi^n/VoKi/2£:o- The motion of the particle in the conservative potential, 
given by the Hamiltonian Hp = p^/2mp -|- hUo |a/o(r)|^, is only governed by the optical 
potential of the pumped mode, while the delayed reaction of the field fluctuations caused 
by the particle also includes the empty modes through their memory operators 



Sn 



/•oo 

/ dre-(''"+^^'^)-/:(r,p;r)/o(r,p;r). (13) 
Jo 

The product of the mode functions /*/o is evolved back in time in the same manner 
as in ( ITTl) . The simple geometry of a particle in an empty Fabry-Perot cavity with a 
transverse pump laser beam is depicted in Figure [H In this case the pumped mode /o 
could be a running wave in y-direction, while the cavity modes would be one or more 
standing wave modes in the x-direction. 

We will discuss in the following section that the cavity memory effect results in 
friction and diffusion effects which add up with increasing mode number. 
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3. The master equation in phase space 

In order to study force and diffusion effects, that arise from the delayed cavity feedback 
on the motion of the particle, it is expedient to translate the underlying master equation 
into its Wigner-Weyl phase space representation [43l HH |45]. Here, the equivalent of 
the motional quantum state p of a particle is given by the Wigner function w{r,p) 
of the phase space coordinates r,p. Position and momentum operators acting on 
the state translate to multiplications and vectorial derivatives by the respective phase 
space coordinates, accordingly. The details of the translation procedure are given in 



Appendix B We treat the particle motion as one- dimensional along the axis of the 
cooling cavity, i.e. we will resort to scalar phase space coordinates {x,p). A realistic 
resonator may often be regarded as one-dimensional if the particle motion perpendicular 
to the cavity axis is sufficiently well described by a classical trajectory, where the 
influence of the cavity is negligible. One may then replace the respective position and 
momentum operators in the master equation by the classical values at time t. This yields 
a parametric time dependence of the mode function, which may be averaged over. A fast 
directed beam of hot particles crossing the cavity, for example, would experience such 
an average feedback interaction during the passage time through the cavity. Particles 
which are trapped in a region overlapping the cavity mode volume, on the other hand, 
could benefit from a larger average interaction time. 

We first translate the master equation ([9]) of a particle in a single pumped mode 



in one dimension. According to the rules given in [Appendix B[ we arrive at a partial 
differential equation for the particle Wigner function with an infinite sum of derivatives 
in X and p. A semiclassical approximation that disregards all derivatives higher than 
second order leads to a Fokker-Planck-type equation 

dtw{x,p) = -d^ {g^{x,p)w{x,p)) - dp {gp{x,p)w{x,p)) + {Dr,^{x,p)w{x,p)) 

+ ^dp {Dpp{x,p)w{x,p)) + d^dp{D^p{x,p)w{x,p)) (14) 

which contains the leading order drift term g^, the force gp, and the diffusion terms D^x, 
Dxpi Dpp acting on the particle. We can therefore easily identify the parameters that 
lead to dissipation and slowing of the polarizable particle. The equation can be split 
into two parts, 

dtw{x,p) = [dtw{x,p)]^^^ + [dtw{x,p)]^.^ . (15) 

Here, the first part is associated with the coherent evolution under the Hamiltonian 
Hp in (|9]). The second part describes the dissipative cavity reaction. The Hamiltonian 
evolution corresponds to the quantum-Liouville-Equation in phase space, as discussed 
in Appendix B It reduces to the classical Liouville equation in the semiclassical limit, 
that is, it evolves the Wigner function according to the classical equations of motion 
under the infiuence of the conservative potential force 

g^^^'^\x,p) = -hUo\a\'dx\fix)\\ (16) 
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It does not contribute to diffusion. We will therefore omit any further discussion of the 
coherent part and focus entirely on the dissipative evolution. The crucial step here is 
to translate the delayed mode function ( ITTi) into phase space. By construction, its time 
derivative fulfills the same quantum-Liouville equation, which we again approximate 
by the classical Liouville equation. The memory operator ( JTOl) then translates into the 
memory integral 



poo 

/ dre-(«+^^)-|/(a:!i,(x,p)) 

^0 



(17) 



where the past position x^^!^(x,p) at time — r is obtained by evolving the phase space 
coordinates (x,p) back in time along the classical trajectory that is governed by the 
conservative force (|T6l) . This guarantees that the delayed mode function expression 
fulfills the Liouville equation. 

In order to assess the friction and diffusion behaviour analytically, we approximate 
the past trajectory in the memory integral by the free motion x^^(x,p) ^ x — pr/mp. 
This corresponds to the above mentioned shearing approximation of Equation ( ITTj) . We 
are allowed to do so, if the average momentum change from the conservative force (fT6ll 
during the memory time I/k is small compared to the momentum p of the particle, that 
is, 

Uoa 



kv 



K 



< 1. 



Here, we used the fact that the derivative of the mode function in (fT6|) gives a factor 
of the order of the wave number, dxf{x) ~ kf{x). We denote the Doppler shift of 
the moving particle by kv = kp/mp. The recoil frequency Ur = hk'^/2mp is the kinetic 
energy transferred to the particle by a single scattered photon, which is a small quantity 
for the large masses considered in this article. In total, the above condition is violated 
when either the particle is so slow or the pump amplitude a is so strong that the optical 
potential of the pumped mode starts to act as a trap for the particle. We do not consider 
this case here, and we approximate the memory integral from here on as 

2 



G{x,p) 



f [x 

mp 



(19) 



Putting everything together, we arrive at the force and diffusion terms of the 
dissipative evolution 



hU^\a\^^t [2tG{x,p)d^ \f{x) 
-hid,G{x,p))dl\f{x)\'] 



D^^'\x,p) = 2feU^\a\''m [dMx,p)]d,\fix)\'' 
D^i'\x,p) = - ffU^\a\'m [d,G{^,p)] \f{x) 



(20) 
(21) 
(22) 



There is no contribution to the coefficients and Dxx in the Fokker-Planck equation 
f[T^ . The second line in the force term fl2U]) represents the second order correction in h 
to the dissipative force, while the first line is of order O (h^)- Again, we can generalize 
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the procedure to N particles in M cavity modes to obtain the dissipative terms as given 
in 



Appendix D 



The Fokker-Planck coefficients ( l20l) and ( |2T1) are the basis for studying the motional 
damping of the particle due to the delayed reaction of the cavity. A closer look at the 
velocity dependence of the dissipative force (l20l) reveals whether the particle is slowed or 
accelerated by the cavity. At low particle velocities v = p/mp the force can be expanded 
to ffist order, 

= gl^'^^x, 0) + P{x)p + O {p') (23) 

where the second term represents the friction force acting on the particle. Motional 
damping is indicated by a negative sign of its friction coefficient P{x) = dpg^^^\x,p — >■ 
0), which is 

If f3{x) < the coefficient can be seen as a momentum damping rate in units of Hz. 
The derivative of the mode function f{x) gives the factor dxf{x) ~ kf{x). This leads 
to the criterion \kv\ < k for the applicability of the low- velocity expansion (!23|) : The 
dissipative force can no longer be described in terms of the friction coefficient (IMl) if the 
particle travels more than the distance 1/k within the cavity reaction and decay time 

The ffist and the second line in flM|) correspond to the friction in ffist and second 
order of the recoil frequency Ur, respectively. 

It is interesting to note that the sign of both contributions to the friction is fixed by 
the cavity-pump detuning A, irrespective of the particular mode function. In particular, 
red-detuned pump frequencies, i.e. A > 0, yield a momentum damping in ffist order 
of the recoil frequency and our result for the friction then agrees with the results in 
[351116], which were obtained in a purely semiclassical approach with Langevin equations 
of motion. 

This represents the intuitive picture of enhanced Doppler-shifted backward- 
scattering of a pump photon, transferring kinetic energy of the order of hwr to a blue- 
detuned and rapidly damped cavity photon. The effect is strongest at the steepest slope 
of and it vanishes at the extrema of the optical potential. Consequently, the 

friction force is significantly suppressed when the optical potential is deep enough to 
trap the particle at a potential minimum. 

The optimal detuning maximizing the magnitude of the ffist order term in is 
A = k/\/3. Putting everything together, the position-averaged friction then scales as 

Y~ 

We find that the momentum damping is a second-order effect in the weak-coupling 
parameter Uo\a\/K,. Note that a very weak polarizability contained in Uo can be largely 
compensated by a sufficiently strong field a. For too strong a the condition of weak 
coupling is not fullfilled anymore which, however, does not necessarily mean that no 



/5 



3^3 Una ^ /u;2 



Ur + 0[^]. (25) 
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cooling occurs. The friction is in addition proportional to the recoil frequency, which 
decreases with growing mass of the particle. In order to achieve a significant slowing 
of the particle it has to interact with the pumped mode over a time of the order 1/13. 
A similar calculation shows that the position-averaged momentum diffusion coefficient 
scales as 

3hmp 

Dpp — 

at low particle velocity. 

The discussion of the friction term in fl24p . that is of second order in u^, is 
omitted here because it is negligible everywhere except for the extrema of the optical 
potential and on resonance A = or |A| ^ /t, where the first order term vanishes. 
At a detuning larger than the cavity linewidth, |A| > k, the second order term anti- 
damps the momentum of the particle. Higher order contributions to the force could be 
obtained beyond the semiclassical approximation in the phase space translation f|T5l) of 
the master equation. Such contributions are suppressed by increasing powers of Ur/i^- 
This justifies, a posteriori, the semiclassical approximation for sufficiently large particles, 
where Ur/ k <^ 1. 

Realistic friction rates can be estimated from the characteristic coupling parameters 
given in Table [1] for a number of large polarizable particles. We consider a Gaussian 
mode with a typical mode volume V = 0.1 mm^ and a damping rate k = 1 MHz, that is 
pumped by an IR laser of wavelength A = 1.5 ^m. A laser with a power of Pl = 0.25 W 
would pump this mode to a stationary photon number |ap = PL/^Khup ~ 10^^. At such 
a power the single-particle weak coupling condition \UQa\ <C n would still be fulfilled also 
for the heavy particle examples in the table. The average friction expression fl25l) thus 
holds and yields the velocity damping rates /3 = 38 Hz, 3 Hz, and 0.8 Hz for the lithium, 
the silica, and the gold cluster in the table, respectively. Although the polarizabilities 
of the three cluster types are roughly the same, the lithium cluster profits here from its 
comparably small mass. Compared to the metal clusters, both the Cgo fullerene and a 
helium nanodroplet of 1000 atoms have a much smaller damping rate of about 100 mHz, 
whereas a single lithium atom reaches 800 mHz because of its light weight. Note that 
the laser and cavity parameters used here could lead to a violation of the weak coupling 
condition by the simultaneous presence of already a few of the heavy clusters in the 
table or of about 1000 fullerene particles, for instance. 

If the weak coupling condition is violated, the friction and diffusion terms derived 
here cease to be valid. However, complementary semiclassical approaches to this problem 
and numerical simulations in the strong coupling regime [35| B6] suggest that the friction 
should generally scale further with the coupling strength beyond the weak coupling 
limit. In general, the effective strong coupling regime \UQa\ > k can be achieved in 
the experiment with large, and thus highly polarizable, atomic clusters in connection 
with a strong pump laser power. On the other hand, the counteracting decrease of 
the recoil frequency with the mass inffuences both the conservative and the dissipative 
forces acting on the particle, as well. 



Uoa 



KUJr 



(26) 
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There are additional momentum diffusion effects due to the absorption and free- 
space scattering of the laser light to be considered, which accompany the diffusion caused 
by the cavity alone. If they are comparable to or larger than the cavity diffusion effect 
they will significantly influence the cooling limit of an ensemble of particles, i.e. the final 
kinetic energy reached by the friction process. The additional diffusion effects will be 
discussed in Section |H A simple damping- diffusion time evolution model for the velocity 
[171 135] may then be applied to estimate the final kinetic energy where the friction effect 
and the momentum diffusion due to the cavity compensate each other. 



\2mp/^ 2mp^{x) ~i\ A, 

Here, we used the fact that the momentum diffusion coefficient (12T|) exhibits the same 
functional dependence of the mode function as the first order friction coefficient in (12^ 
in the limit of small particle velocity. The expression ( 127]) depends on the detuning and 
ceases to be valid in the extreme cases A = 0, A ^ k, where the first order friction 
vanishes. In case of the above optimal detuning the final energy would be hK/2\/^. 
This recovers the result of the single-particle cavity cooling limit and agrees with the 
results of wavefunction simulations in [10]. For an experiment it suggests the use of 
very good cavities to obtain faster cooling by larger field amplitudes together with lower 
temperatures k. 



Table 1. Light coupling parameter table for different polarizable particles and IR 
laser light with a wavelength A = 1.5 /im. The polarizability x and the absorption cross 
section a a at this wavelength are listed for particles of different mass m, with a focus on 
heavy-sized clusters of gold, lithium and fused-silica. Their optical properties are taken 
from [48l|4^. For the lithium atom and the superfluid Heiooo droplet the polarizability 
is estimated by its static value per atom [50l [51] . Both He and Li are transparent to IR 
light and we omit their absorption here. For the Cgo-fuUerene we extract both x ^^^d 
a a from the spectra in [52j . The derived values of the recoil frequency ujr , the coupling 
constant Uq, the photon absorption rate 2ja, and the Rayleigh scattering rate 2^s are 
given in units of MHz, which corresponds to typical, experimentally achievable cavity 
damping rates k. The cavity mode volume is set to ^ = 0.1 mm'^. 
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• 10-6 


3.3 
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4.8 • 10-9 


2.7 • 10-13 



3.1. Fabry- Perot cavity 

For the instructive example of a one-dimensional Fabry-Perot cavity with mode function 
f{x) = cos kx the memory integral (1191) and thus the dissipative force and diffusion terms 
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can be calculated explicitly for any given velocity. Focussing on their dependence of the 
particle velocity, we take the average over the position coordinate and obtain 

^(^P.) = - 2»/^^^ (28) 

with the complex parameter = k + iA. The force term is given in leading order of 
the recoil frequency Ur- Both terms scale identically with the effective coupling to the 
cavity field ?7o|q^|- Moreover, they vanish in the Doppler resolution limit kv <^ k,A, 
where the delayed cavity reaction is too slow to affect the particle's motion. The velocity 
dependence of the position-averaged dissipative force is plotted in Fig. |2]for an effective 
coupling \UQa\ = 0.1k and different detunings A. One clearly sees the linear frictional 
behaviour of the force term in a small window around zero velocity, where the detuning 
gives the optimal slope. On the other hand, the range of velocities that are 
damped by the dissipative force term increases with growing detuning A. In the limit 
of large Doppler shifts \kv\ ^ n the damping effect is maximized for A ^ \kv\. Fast 
particles are therefore best addressed with a detuning far off the cavity resonance. Note 
that in order to combine large capture ranges and strong friction for low velocities one 
could think of a time modulation of the detuning or the simultaneous use of various 
cavity modes operated at different detunings. 

Note that our approximation of the memory integral (fT9|) formally fails in close 
vicinity of the point kv = 0, where the condition (fTSl) is violated. Although, strictly 
speaking, the discussed friction force is therefore not valid for such small velocities, we 
expect the error to be negligible for typical pump amplitudes that do not by themselves 
already create an optical potential that is strong enough to dominate the particle motion. 

We also want to point out that the positivity of the average momentum diffusion 
( l29l) does not guarantee the positivity of the diffusion coefficient at every point {x,p) in 
phase space. This is in fact not the case in our semiclassical approximation, as we will see 
below. The positive semidefiniteness of the diffusion matrix (-Djj) is required to interpret 
the semiclassical Wigner function evolution equation that is derived from our master 
equation as a classical Fokker-Planck-Equation [53]. Only then does it make sense to 
translate the phase space equation into an equivalent classical Langevin equation for 
the position and momentum coordinates. However, in the Fabry-Perot case considered 
here, the momentum diffusion coefficient fl2T]) reads 



+ [2kvY\ 

- 2kv (3?e + {2kvY) sm{2kx) cos(2A;a;)] (30) 

As demonstrated in Fig. |3] it is not positive everywhere for nonzero velocities. The 
diffusion coefficient takes on small negative values close to the minima of the optical 
potential, but it is still dominantly positive for slowly moving particles, \kv\ < k. In the 
case of faster particles, where the negative and positive regions almost compensate each 
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Figure 2. Position-averaged dissipation force (pS)) plotted versus the particle velocity 
for different cavity-pump detunings, using an effective coupling strength \Uoa\ — O.Ik. 
The velocity is expressed in terms of the Doppler shift kv in units of the mode damping 
K, and the force is plotted in units of the photon momentum transfer rate {hk)K. The 
solid line represents the detuning A = k/ -s/S, which gives the the optimal friction 
coefficient (p5)) . that is, the steepest slope around kv = 0. The dashed and the dotted 
line represent the higher detuning A = 2k and A = 5k, respectively. Although they 
exhbit a weaker damping of small velocities, their damping range of velocities grows. 




other, the momentum diffusion effect vanishes rapidly. We attribute the negativity to the 
non-Markovian nature of the underlying master equation iQ. The memory effect of the 
cavity, and in particular the accompanied correlation between the particle motion and 
the cavity field fluctuations [54j, impede a completely positive Lindblad form [55 | [56 | 157] . 
The effect is however small and negligible in the practical considerations on the motional 
damping of large polarizable particles discussed in the following. 



3.2. Confocal cavity 

As mentioned in Section [21 the effective friction force should increase in the presence 
of M near degenerate empty cavity modes /i, . . . , /a/ in addition to the pumped mode 
/o- Again, this can be best studied after translating the corresponding master equation 



(1121) into phase space; this is a special case of the general formalism in [Appendix D 
Disregarding the Hamiltonian evolution and considering only one dimension of motion 
the semiclassical dissipative force and diffusion coefficients read as 

M 

g^'''''\x,p) = \Uma\^^t [2zGUx,p)dJ*{x)fUx) 



m=0 

n2 f* 



-h{d,GUx,p))d'J*ix)fUx)\. (31) 



M 



Dlf\x,p) = 2ff \Uma\'^c [d,GUx,p)dJ*{x)fUx)] , (32) 



m=0 
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Figure 3. Momentum diffusion coefficient ([50]) plotted as a function of position 
for several fixed velocities. The particle's position x is expressed in terms of the 
dimensionless parameter kx in the cavity field, and the diffusion coefficient is given 
in units of {hk)'^K,. The curves correspond to the momentum diffusion expressions for 
the fixed Doppler shifts kv = (dotted line), kv — 0.5k (dashed line), kv = k (dash- 
dotted line), and kv = 5k (solid line). The negative part of the diffusion coefficient 
grows with increasing modulus of the velocity, until it practically compensates the 
positive part. Negative velocities give the same results due to symmetry. Only the 
diffusion for kv — is non-negative everywhere. 



M 

Df;'\x,v) = -^'Y. \Uma\^^t [dpGUx,p)dJ^{x)frn{x)] , (33) 

where we introduced the memory integral for each mode in the same free-motion- 
approximation as before 

Gm{x, p) = r dr e-(-'"+^'^'") (x - ^) f, (x - . (34) 

Both the force and the diffusion coefficient are sums over terms similar to the single- 
mode case (120]) - (122]) . The only difference is that the modulus square of a single-mode 
function is replaced by the cross-products f^fm and the complex conjugates. In the 
limit of small velocities the friction coefficient now reads 

M ,jj |2 

= - E ^ifrT. V i^rn \d. (/o(x)/.(x))r 

+ l^ ''\~^} \dlU;{^)U^))\') ■ (35) 
2mp + A^^ ' ' y 

Again, the sign of the low-velocity friction is independent of the mode functions. It is 
solely determined by the signs of the mode detunings with respect to the pump 
frequency. 

The multi-mode friction force can be illustrated with the practical example of a 
confocal standing wave cavity. In order to apply our one-dimensional description here, 
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we neglect all dynamical effects due to the transverse mode profile, and we treat the 
transverse coordinates z as parameters to arrive at effectively one-dimensional modes 
along the cavity axis x. In the end, we take the average of the resulting friction coefficient 
( 135|) over the off- axis coordinates y,z. By setting the averaging area much larger than 
the cavity mode extensions, this approximation is relevant for experiments where a 
directed molecular beam crosses the cavity to obtain transverse motional cooling. The 
particle ensemble is assumed to be sufficiently dilute so that only one particle passes the 
cavity at any given time. In a typical case, when the mode waists are much larger than 
the laser wavelength, the transverse dissipation and dipole forces in z can be safely 
neglected. 

We consider a confocal cavity with two circularly symmetric mirrors of sufficiently 
large aperture, where the mode functions are well described by Laguerre-Gaussian 
standing- wave modes [CT 136] . They are characterized by a longitudinal mode index 
n e N and two transverse mode indices m, £ G Nq that yield the resonance wavenumbers 



The distance d of the mirrors equals their radii of curvature. Hence, for a given 
fundamental mode (no,0,0) there exists an infinite number of degenerate higher-order 
transverse modes (n < nQ,m,i) at the same resonance frequency. In practice, the 
extended transverse profiles of higher order modes, however, limit the degeneracy as the 
mode damping K,n,m,i increases with growing m, i. This can be expressed in terms of the 
effective mode waist [36] 



It grows with respect to the waist wq of the corresponding fundamental mode (no, 0, 0). 
The maximally achievable waist without a significant loss of cavity finesse then serves as 
a rough estimate of the number M of degenerate modes: If a confocal cavity for example 
tolerates an effective waist without considerable losses, that is a times larger than the 
waist of the pumped fundamental mode, then the transverse mode index expression 
2m + i + 1 may have values up to a^, where both m and i can vary independently. 
This results in a degeneracy of the order of modes which contribute to the average 
friction. These modes include both cosine and sine standing wave modes for n even 
and odd, respectively. The advantage of this is demonstrated in Fig. HJ where the one- 
dimensional friction coefficient ( l35l) is plotted as a function of x-position around the 
center of a confocal standing-wave resonator for different contributing mode numbers. 
We consider the fundamental mode no = 2 ■ 10^ of the confocal cavity of size d = 10 mm 
to be coupled to a Gaussian running wave mode /q in y-direction with the same waist 
Wo, following the scheme in Fig. [H The y-mode is strongly pumped by an IR laser 
detuned with respect to the cavity to A = n/y/S, and the mode wavelength is given 
by A = 2d/{nQ + 1/2) = 1 /xm. We assume that all degenerate modes share the same 
common linewidth n. The pump amplitude a is fixed to a value that corresponds to the 
effective couphng UnoOi = O.Ik between the pumped Gaussian mode and the fundamental 




(36) 



Wn,m,i = Wo\/2m + i +1. 



(37) 
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Figure 4. Plot of the modulus of the multi-mode friction coefficient ((35)) of a confocal 
resonator versus the position kx along the cavity axis. The transverse coordinates y 
and z are averaged over a large area with respect to the waist of the contributing cavity 
modes. The pump mode is given by a Gaussian running wave mode in y-direction. See 
the text for numerical details. The dotted curve represents the friction related to the 
coupling of the pump mode to a single fundamental mode of the cavity only, whereas 
the pump mode couples to M = 15 and M — 55 degenerate modes (including the 
fundamental mode) in the dashed and in the solid lined plot, respectively. The friction 
coefficient is negative in all cases, and it is plotted in units of 10^^ times the linewidth 
of the cavity modes k. 



cavity mode. We set the recoil frequency of the particle to u.,. = 0.001k in this case, 
which corresponds to a mass mp = 1000 amu and a cavity linewidth k = 1.25 MHz. The 
x-derivative of the running-wave mode /o may be neglected in the multi-mode friction 
formula (|35l) . i.e. the friction is only caused by the coherent scattering of pump photons 
into the degenerate confocal cavity modes. The transverse coordinates y and z are each 
averaged out numerically over a distance of 4 times the fundamental waist off the central 
cavity axis. Fig. H] shows that the average friction not only increases with growing mode 
number, as expected, but it also does not oscillate to zero anymore. The latter is due 
to the equal presence of both sine and cosine standing wave modes. 

In practice, all degenerate high-finesse modes, which are supported by the confocal 
resonator, contribute to the friction. The dependence of the friction coefficient fl35|) 
on the supported degeneracy is depicted in Fig. [5l Here, the one-dimensional friction 
of Fig. m is in addition averaged over x and plotted as a function of the number of 
contributing degenerate modes. The mean friction experienced by each member of 
the particle ensemble crossing the cavity increases with each higher-order mode that 
is supported by the cavity. However, Fig. [5] does not indicate a linear growth of the 
average friction with the degeneracy. This can be explained with the multi-mode friction 
formula (l35|l . Each empty cavity mode fm enters there in a product with the Gaussian 
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Figure 5. Modulus of the position-averaged friction coefficient of Fig. |4] plotted as a 
function of the number M of contributing degenerate modes in the confocal cavity. The 
first data point represents the friction related to the coupling of the Gaussian pump 
mode to the single fundamental mode [tiq = 2 • 10"*, 0,0) of the cavity. It corresponds 
to the average of the dotted curve in Fig. 01 All data points in this plot are given 
relative to this value. The rightmost data point contains the friction contribution of 
M = 406 degenerate modes including the fundamental one. 



pump mode /o, which has a fixed spatial extension given by its waist wq. Then again, 
the transverse profile of higher order modes stretches over a larger area off the cavity 
their effective waist ( 1371) grows. Consequently, their overlap with the pump 
mode descreases, and the multi-mode friction coefficient saturates for a sufficiently large 
degeneracy. 

4. Light absorption and scattering 

So far we have assumed the particle to have a real, scalar and positive linear response 
X to the optical cavity field, i.e. to be an ideal high-field-seeker with no photonic loss 
channels such as absorption or scattering into free-space. This has led to the derivation 
of an effective master equation in Section [2] and yielded the analytical expressions for 
friction and diffusion in Section [3l In practice this condition holds fairly well for two- 
level-atoms in optical fields far red-detuned from their internal resonance [58j. It is also 
a good description of large molecules or clusters with sufficiently low absorption, where 
the large number of internal degrees of freedom in general does not allow us to formulate 
a simple few-state coupling model. Nevertheless, many theoretical approaches towards 
dissipative cavity cooling of molecules [HI [59l [60] start from an effective two-level- 
description of the molecule's internal dynamics, which is then adiabatically eliminated 
in the off-resonant limit to arrive at a real polarizability again. In such a formulation one 
can easily incorporate spontaneous emission of photons into the vacuum, which yields 



Master Equation for the Motion of a Polarizahle Particle in a Multimode Cavity 19 



a free-space Rayleigh scattering rate 7^ in the off- resonant limit. 

In the following we discuss the influence of scattering and absorption at large 
particles in terms of experimentally accessible scattering and absorption cross sections. 
We introduce additional Lindblad terms in our master equations to account for these 
effects, and we translate them into phase space in order to reveal their influence on 
dissipation and momentum diffusion. The dominant effect is caused by the pump 
photons and we neglect the influence of field fluctuations or thermal photons. 



The derivation of the additional Lindblad terms is given in Appendix C where 
we apply the same weak coupling model as used in Sec. |2] to couple the pump field to 
internal degrees of freedom of the particle (to describe absorption) and to free-space 
vacuum modes of light (to describe Rayleigh scattering). The weakness of the coupling 
in these cases is a reasonable assumption for the large particles and the off-resonant 
pump light considered here. This is illustrated by the examples in Table [H where the 
the single-photon absorption and scattering rates 7a, 7s are already significantly smaller 
than the single-photon coupling Uq to the cavity mode. 

4.1. Absorption 

In order to achieve dissipation and eventually cooling of a molecular ensemble we have 
to avoid absorption which would lead to heating. Dipole allowed transitions in many 
organic molecules as well as plasmon resonances in metal clusters are typically situated in 
the UV-VIS range. This makes it advisable to pump the cavity with near-infrared light, 
i.e. the extreme tail of the absorption spectrum. Each particle may then be regarded 
as a thermal bath of harmonic oscillators with a linear and approximately constant 



coupling to the external cavity field modes (see Appendix C). When the particles remain 
internally cold during their residence time in the cavity, we are allowed to trace over the 
internal oscillators to obtain the Lindblad term for the pumped cavity mode 

[dtPpU. = 7a|«r (2/(r)ppr(r) - {pp, |/(r)|^}) . (38) 

This term has to be added to the master equations and (|T2l) in order to account 
for the momentum diffusion that is caused by the recoil of absorbed pump photons. 
The absorption rate per photon 7^ can be expressed in terms of the cross section da as 
27a = caa/V. 

4-2. Rayleigh scattering 

The elastic scattering of pump photons into free-space can be accounted for by the 
standard Lindblad term 

Msca =27>r / d'uNiu)Kppkl-^s\a\' {pp,\fi'r)\'} (39) 

J\u\=l 

K = e-'^^^-'/'f{r) (40) 



As shown in Appendix C, this expression is obtained by modelling the space outside 



the cavity as a band of vacuum modes with arbitrary detuning in a master equation of 
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the form fll2p . The couphng parameters are then replaced by the Rayleigh scattering 
rate 7^ and the angular scattering distribution N{u). Any Doppler shift related to 
the scattering by moving particles can be neglected since we are only interested in the 
additional momentum diffusion that is caused by the random recoil of the scattered 
photons. In a one dimensional situation, where only the particle motion along the 
cavity axis x is taken into account, one can resort to an effective mode function f{x). 
In this case the scattering Lindblad operators are replaced by 

\,^ = e~^-^-^-/j{x), (41) 

with the component of the vector u along the cavity axis. Recall our sans-serif 
notation for the position operator x that distinguishes it from the scalar value x. We 
obtain the scattering rate 7^ of pump photons from the Rayleigh scattering cross-section, 
7^ = cas/V. 

Both rate constants 7a and 7^ depend on the dielectric properties of the particle. 
In case of atomic clusters that are sufficiently large to be modeled as dielectric spheres 
- but still small compared to the light wavelength - the rates can be extracted from the 
complex polarizability. It reads as [611 [62] 



with R the cluster radius and e{X) its complex dielectric function at the wavelength A 
of the applied light field. This results in the absorption cross section 

2% ^ 27iR. fe(X)-l\ ,^ , 

section 
e{\) - 1 



and in the Rayleigh scattering cross section 

2 



Svr 2TiR , ^ 



.(A) + 2 • ^''^ 

Knowing the dielectric function e(A) at the laser wavelength one can compute the 
scattering and absorption rates, as has been done for the exemplary cluster particles in 
Table [H Although the resulting rate values there are by orders of magnitude smaller 
than the cavity decay rate, the momentum diffusion that is associated to the absorption 
and the Rayleigh scattering of pump photons will generally surpass the diffusion related 
to the cavity modes, as we will discuss below. In the particular case of sub-wavelength 
metal clusters the absorption rate is much larger than the elastic light scattering rate, 
7a ^ Is- Consequently, the absorption effect will dominate the momentum diffusion 
along the direction of the pump laser, and the scattering effect may be neglected. This is 
however not the case in a two-dimensional setting like the one in Fig. [H where the pump 
mode and the empty cavity modes are orthogonal to each other. Here, the absorption 
effect merely causes momentum diffusion along the pump direction while the diffusion 
along the x-direction is mainly related to the Rayleigh scattering of photons out of 
the pump beam. Such an experiment would therefore benefit from the much smaller 
scattering cross sections. 
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4-3. Realistic dissipation rates 

We have now properly included light absorption and scattering into our master equation 
model which, after translating those terms into phase space, give us a realistic description 
of the friction and diffusion dynamics of large particles in a cavity. The translation 
mainly leads to an additional momentum diffusion that counteracts the envisaged 
frictional momentum damping in dependence of the respective absorption and scattering 
rates. 

A Rayleigh scattering event can be divided into the momentum recoil due to the 
absorption of a photon from the pump mode, and the recoil caused by the reemission of 
the photon. The former part is identical to the effect of genuine photon absorption, 
i.e. the scattering rate 7^ can be added directly to the absorption rate 7^ in the 
absorption Lindblad term ( l38|) . In a one- dimensional situation this yields the non- 
negative momentum diffusion 

D(f = h\^, + 7.) \adj{x)\^ (45) 

along the cavity axis x. Only x-derivatives of the mode function are involved. 
Consequently, this diffusion term can be suppressed for instance by using a large waist 
of the pump mode perpendicular to the axis x of a cavity with a number of empty 
modes. The latter will cause the desired friction in x according to (135|) . while the pump 
mode fulfills d^f ~ and thus minimizes the additional diffusion term ( 145|) . In case of a 
running wave pump mode f{x) ^ M directed along the x-axis the absorption also leads 
to a net force 

gf'''\x,p) = 2M7a + 7s)|«r3m [f{x)dj*{x)] (46) 

that acts on the particle in the preferred direction of the photon flux. It vanishes for a 
standing-wave pump where the same amount of photons hits the particle from one side 
and the other. 

The reemission part of the scattering results in the diffusion coefficient 

D(7)(x,p) = [^{ul)\f{x)\' + 2{u.,)'^m [f{x)dj*{x)]) . (47) 

The angular scattering distribution N{u) enters through the average x-component of 
the scattered photon momentum {ux) = f d'^u N{u)ux, and the averaged square (u^). 
The second term in the diffusion coefficient (l47j) is caused by a combined anisotropy of 
the photon flux of the mode, /(x) ^ M, and of the scattering, (u^) 7^ 0. The latter also 
implies a net force 

^?f^)(x,p) = 2%,5^(n.)|a/(x)r (48) 

acting on the particle. Both anisotropic forces f H6|) and fHHj) are independent of the 
particle velocity and therefore do not contribute to friction. 

If the particle scatters photons isotropically into free space, the average (u^) = 
and {uD = 1/3. It is possible to adjust (m^), i.e. the diffusion due to Rayleigh scattering, 
by using linear polarized pump light from a mode perpendicular to the cavity axis x, for 
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instance. When the polarization axis is parallel to the cavity axis, one obtains {u^) = 
and (u^) = 1/4 from the corresponding dipolar scattering pattern. 

The diffusion terms caused by the absorption and scattering of pump light usually 
govern the total momentum diffusion of the particle as they exceed the momentum 
diffusion effect of the cavity modes discussed in Section |3l In order to compare the 
respective diffusion terms with the diffusion of the cavity, as given by the low-velocity 
term (|26|) we need to average the expressions fH5l) and (H7|) over position. We assume 
isotropic scattering for the latter term. The average diffusion due to absorption then 
scales as (?7o/^)~^(7a + 7s)/^ times the cavity diffusion coefficient, while the scattering- 
induced diffusion scales as (f/o/fi;)~^7s/K times the cavity diffusion coefficient. The 
final kinetic energy that can be reached through the frictional damping process roughly 
increases by the same factors. Regarding the example of heavy clusters listed in Table [T] 
the increase can be dramatic. For the lithium cluster and the gold cluster the diffusion 
caused by absorption surpasses the cavity diffusion by a factor of the order of 10'^. 
However, as already mentioned, the absorption term can be suppressed by using a 
pumped mode that is oriented perpendicularly to the axis of interest x. In this case, 
only the diffusion caused by scattering has to be taken into account, which is of the same 
order of magnitude as the cavity diffusion for all the heavy clusters listed in Tab. [1] Due 
to the high transparency of silica clusters and helium droplets in the IR regime, they 
are not influenced by the absorption-induced diffusion even at high masses. 

5. Conclusions 

In this work we have presented a full quantum treatment in terms of master equations 
of the cavity-assisted motional damping of a polarizable particle. In the limit of 
weak coupling to a strong laser field we could derive closed expressions for both the 
velocity-dependent force and the momentum diffusion related to the delayed reaction 
of high-finesse cavity modes to scattered pump photons. On the single particle level, 
our phenomenological treatment of the particle-light interaction gives insight into the 
challenges and peculiarities of optical cooling experiments with ensembles of mesoscopic 
objects such as large molecules or atomic clusters. Though the weak polarizability of 
those heavy particles can be compensated by a sufficiently strong pump field intensity, 
the friction rate has an upper bound of the same order of magnitude as the recoil 
frequency. This can be fairly slow but our model predicts that it can also be enhanced 
by orders of magnitude when using multimode resonator geometries such as a confocal 
cavity setup. We therefore believe that the presented scheme could prove to be useful in 
practice, in view of its generality and the lack of alternative cooling schemes for freely 
moving massive particles. The large multimode enhancement might even be exploited 
in future optomechanics experiments. For higher temperature sources containing fast 
particles we showed that the capture range and region of optimal cooling can be tuned 
by a suitable choice of detunings and stronger pump. Ultimately, a time- dependent 
detuning can be used to optimize the capture range and the final temperature for a 
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given interaction time. The final temperature is however also sensitive to the absorption 
and the free-space scattering cross-sections of the particle, which should be as small as 
possible. 

Both the frictional part of the force as well as the negativities in the momentum 
diffusion coefficient are a direct consequence of the memory effect of the cavity, i.e. its 
delayed reaction towards the field fluctuations caused by the presence of the particle. 
The memory effect, and in particular the consequential negative diffusion, would play a 
more major role in the dynamics of the particle beyond the weak coupling limit assumed 
here. A strictly semiclassical description in terms of an effective Langevin equation of 
motion for the particle would then have to be replaced by coupled Heisenberg equations 
for the particle and the field. It is, however, still an open problem to obtain an effective 
description in terms of a master equation (without a harmonic approximation of the 
optical potential) for the particle in the strongly coupled regime. 
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Appendix A. General master equation for particles in M modes 

We generalize the derivation of the weak coupling master equation for a particle in a 
pumped cavity mode from Section [2] to A^ particles in M modes. A single laser frequency 
ujp is used for all stationarily pumped modes. We allow for different pump amplitudes 
ctm, detunings A.^ = ujm — i^p, and damping rates Hm- In practice one would single out 
one or at most a few almost degenerate modes out of many to be strongly pumped by a 
single laser, as described by the master equation (fT2l) . It is a special case of the equation 
presented in the following. 

The mode functions fm{.f) determine the spatial structure of the quadrature field 
operator 



which enters the field coupling Hamiltonian ([2]) at each particle's position operator r^. 
For simplicity, we assume the polarizability x to be constant for all modes. In the 
displaced and rotating frame with respect to the pump the master equation of particles 
and modes reads 




M 



(A.l) 



M 



dtp= - T [He + Hp H/, p] + ^ Km {2amp3l, - {a. 




(A.2) 



with the cavity, particle, and interaction Hamiltonians 



M 




(A.3) 



m=l 
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TV / 2 M \ 

= E + E ^ Umnal,aJ*^{r,)frXr,) (A.4) 

j=l \ m,n=l / 

N M 

Hi = U^nQ*m^nf*mi.rj)fn{rj) + h.C. (A.5) 

j=l m,n=l 



The symmetric coupling matrix f/^n = —Xy^m^n/VmVn/'2.So describes the scattering 
of photons at a particle between mode m and n. The dispersive coupling between the 
strong pump fields is responsible for the optical potential in the particle Hamiltonian Hp, 
while the Doppler-shifted scattering of pump photons into a detuned cavity mode at a 
moving particle may cause dissipation of motion. The latter comes from the interaction 
term H/ after adiabatically eliminating the cavity modes. The first weak coupling 
condition in Section [2] stays the same, and has to hold for all modes, \Umn\ ^ i^m 
Vm, n. The second weak coupling criterion however poses a further restriction on the 
particle number, \NUmnOtn\ ^ Wm,n. This is because each particle may scatter 
photons into the same mode m independently, and the weak coupling limit requires the 
cumulative photon fiuctuations of all particles to be smaller than one. 

We may follow the procedure given in Section [2] adapting the Liouvillians there to 
the corresponding terms in f lA.2p . This yields the master equation 

M N 

<9tpp = - - [Hp,pp] - ^ 'YUi.mU.mnia*ian[fifm{rk),gmnPp]+h.C.), (A.6) 

£,m,n=l k=l 

with the A^-particle memory operator 

poo N 

gn.n = / dr e-('^-+^^-)^ f*mfn (r„ p,; t) . (A.7) 

,=1 

The backward-evolved product of the mode functions is given by 

fUn (r, p; r) = e-^^^--/V:.(r)/n(r)e^"--/^ (A.8) 

Appendix B. Phase space translation rules 

The Wigner-Weyl phase space formulation is a way to represent observables and density 
operators of a motional degree of freedom in terms of functions over the phase space 
variables (x,p). The Weyl symbol of an observable A is given by the invertible integral 
transformation 

A{x,p) = j dse'^''^(x-^ k x^^, (B.l) 
and the Wigner function of a state p is its normalized Weyl symbol 

Mx,j9) = ^ /" <\se^^'"^{x-'-\p\x+'-). (B.2) 
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dt + dx — 
mp 



It can be seen as the quantum generalization of the classical phase space density. Given 
the Hamiltonian H = p^/2mp + ^(x) the Wigner function evolves according to the 
quantum-Liouville-Equation |15] 

The nonclassical behaviour entirely comes from the right hand side, and it represents 
diffraction effects at the potential. In case of states p) where the quantum coherence 
is negligible, the semiclassical approximation is applied by dropping all higher-than- 
second-order derivatives of the Wigner function. One obtains a Fokker-Planck-type 
equation f[T^ with the drift terms g^, Qp and the diffusion terms D^x, Dpp, and D^p. In 
a coherent time evolution governed by the conservative force term gp{x,p) = —dxV{x) 
there is no diffusion according to the quantum-Liouville-Equation. This is not the case 
in a general master equation anymore, where the translation into phase space may yield 
a momentum-dependent dissipative force gp{x,p), as well as nonzero diffusion terms. 
The Weyl symbols of the products AB, BA of two operators A, B read 



{AB){x,p) = A[x 



iBA){x,p) 



B{x,p), 



(B.4) 



(B.5) 



where the Weyl symbol with differential operator arguments acting to the right is defined 
via the Taylor expansion 



A 



E 

n,m=0 



J-^A("'«)(x,p) (-^dp 
mini V 2t 



dx 



(B.6) 



The Weyl symbols of operator functions /(x) and g{p) are simply given by f{x) and 
g{p), respectively. 



Appendix C. Lindblad terms of photon absorption and scattering 

In order to derive a Lindblad term that describes the momentum diffusion due to 
absorption or elastic Rayleigh scattering of pump photons we utilize and modify the 
master equation f|T2|) of a particle in one pumped and M empty cavity modes. We 
assume here the same weak coupling assumptions already made in Sec. [2J 

In the case of photon absorption we substitute the empty cavity modes by M ^ 1 
internal undamped oscillators, i.e. we set /m(^) = 1 and = for all those modes. 
This follows the textbook analogy of damping, where an oscillator mode is weakly 
coupled to the infinite free-space heat bath [531 EH E3|- Consequently, we neglect 
any heating effect of the internal bath which would lead to additional diffusion by 
thermal radiation. This is a valid approximation for particles with sufficiently large 
heat capacities, as considered in this paper. We model the internal bath as a continuum 
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of oscillators with arbitrary detunings from the pump and replace the mode summation 
in ( 1T2|) by the integral 

dArf(a;p + A). (C.l) 

m=l 

Here, d{u) denotes the internal density of states. The rotating wave approximation 
around the optical pump frequency Up admits to set the integral bounds to ±00. In 
addition the coupling constant times the density of states is assumed to be flat around 
the pump frequency, 

U{iOp + A)2rf(wp + A)^ U{upfd{ujp) =: 717,. (C.2) 

The constant 27^ is the phenomenological rate of absorption events. Plugging everything 
into Equation ( fT2|) yields the Lindblad term 

[^*Pp]abs = 7a|«P (2/(r)ppr(r) - {pp, \f{r)f}) . (C.3) 

It represents the momentum diffusion caused by the absorption of photons from a single 
pump mode /(r). 

A similar calculation can be performed for Rayleigh scattering of light. Starting 
again from f |T2l) . we replace the empty modes by the undamped free-space plane waves 
/fe(^) = e**''^ with wave vector k = uu/c. We disregard the Doppler shift of the 
elastically scattered pump photons, since we focus on the diffusion effect. The memory 
operator f|T3l) then becomes 

POO 

/ dTe-'^^f{r)e~'''-' ^n5{A)f{r)e-''^^''-'/'. (C.4) 
Jo 

The principal value due to the semi-infinite integral is dropped as it would only lead to 
a negligible Lamb shift term [65l |66]. Given a phenomenological scattering rate 7^ and 
the distribution of scattering angles N[u), one arrives at the Lindblad term (l39l) . In 
the case of Rayleigh scattering of linear polarized pump light the distribution is given 
by a dipole scattering pattern N{u) oc sin 6 with 6 the polar angle with respect to the 
polarization axis. 

Appendix D. General friction and diffusion terms 

We translate the dissipative part of the general master equation for particles in 



gk 



M modes from Appendix A into phase space, using the translation rules derived in 
Appendix B We restrict the calculation to the one-dimensional case along the cavity 
axis X. Note that the phase space of A^ particles is 2 A^- dimensional then, with 
coordinates {xi,pi, . . . ,xiy,piy). A tedious but straightforward calculation leads to the 
following force and diffusion terms in the semiclassical limit: 

The force consists of a first order and a second order contribution in h, 
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Both contributions are expressed in terms of the memory integral per particle 

GUx,p) = rdre-('^-+^^-)V; (x - ^) fn (x~—). (D.2) 
Jo V ^pJ V ^pJ 

The past position in the argument of the mode functions is approximated by the free 

trajectory x + pr/mp, as it is done in Section^ 

The dissipative force terms now read 



M 

gpJ{Xi,Pi, . . .,Xn,Pn) =2h ^ UlmUmn^t [ia^an {dx^fKXk)fm{Xk)) 

£,m,n=l 
N 

^ ] Crmn{Xj,Pj) 



N 

X 

i=i 

M 



(D.3) 



9pJiXl,Pl,---,XN,PN) = -fi^ UimUmn^^^ [<y*eOin{dlji{Xk)fm{Xk)) 

£,m,n=l 

X dp^GmniXk,Pk)] ■ (D-4) 

In the first order term all particles add up homologously via their memory integral to 
the force acting on particle k. The second order force is however only related to the 
field fluctuations caused by the particle it acts upon. Finally, the nonzero parts of the 
diffusion matrices read 

M 

DpkPjixi,pi, . . .,xn,Pn) = 2/i^ ^ UemUmn^t [a|a„ {dxje{xk)fmixk)) 

£,m,n=l 

X dx^Gmn{Xj,Pj)] , (D.5) 
M 

Pl, . . . ,Xn,Pn) = - ^ UimUmn^t [a}an{dxJiiXk)fmiXk)) 

£,m,n=l 

X dpPrnn{Xj,Pj)\ . (D.6) 

In both diffusion and first order force the particles couple to the field fluctuations caused 

by the other particles in the same way as they couple to their own field fluctuations. 

This indicates a strongly correlated system, where one naturally expects collective effects 

beyond our weak coupling limit, such as self-organizing phase transitions [T8l [13] . 
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